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Qs I models of electroweak baryogenesis which produce a highly inhomogeneous distribution of 
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can be ruled out whenever a fraction /^3 x 10 ^(100 GeV/T)^ of nucleated bubbles of 
broken-symmetry phase contributes ^10% of the baryon number within a horizon volume. 
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In this letter we discuss how the sensitivity of big bang nucleosynthesis to the baryon- 
to-photon number (rj) and its spatial distribution could be utilized to probe electroweak 
physics in a new manner. Models of electroweak baryogenesis have so far concentrated on 
producing the presently observed value for r]. In this paper, we point out a new constraint 
on models that result in a highly inhomogeneous distribution of r]. We will show that 
micro-physical processes that generate fluctuations in ry on sub-horizon scales for epochs 
corresponding to temperatures T ^ 1.5 TeV may be subject to nucleosynthesis constraints. 

A long-standing problem in astrophysics is the explanation for the apparent baryon 
number asymmetry in the universe. Reference [1] provides an overview of this problem 
and the attempts to solve it. Ever since the work of Sakharov [2] an explanation has been 
sought for the baryon number asymmetry in C, CP, and baryon- number violating processes 
in environments associated with departures from thermal and chemical equilibrium in the 
early universe. However, any net baryon number generated at very early epochs in the his- 
tory of the universe (e.g., via C and CP violating, nonequilibrium baryon number violating 
decay of heavy X and Y bosons associated with Grand Unification) will probably, though 
not necessarily inevitably, be erased by subsequent anomalous electroweak processes [3]. 
Regeneration of baryon number could then occur during a first order cosmic electroweak 
symmetry-breaking phase transition [3] . It is not clear if adequate baryogenesis could be 
achieved with a minimal Weinberg- Salam model without implying a Higgs mass below 
the present experimental lower bound (see, however, reference [4]). Several plausible ex- 
tensions of the minimal standard model, such as multi-Higgs models or supersymmetric 
models, could lead to significant baryon number generation at this epoch [5,6,7]. A re- 
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view of baryogenesis associated with first-order electroweak phase transitions is given in 
reference [8]. 

Since a temperature dependent nucleation rate is a generic feature of first order phase 
transitions, we expect some supercoohng in a primordial electroweak transition and the 
concomitant generation of distinct bubbles of low temperature phase. These bubbles of 
broken phase grow until they coalesce. As the bubble walls propagate toward coalescence 
the universe is out of thermal and chemical equilibrium in the vicinity of the walls. These 
nonequilibrium conditions, together with baryon number violating anomalous electroweak 
interactions and C and CP-violation, provide all the necessary ingredients for baryogenesis. 
The necessity of this baryogenesis occurring in the inhomogeneous environments engen- 
dered by bubble nucleation and coalescence ultimately may lead to an inhomogeneous 
distribution of ry and, hence, entropy-per-baryon. 

In most of these baryogenesis scenarios the final distribution in ry is probably too 
homogeneous to affect nucleosynthesis. However, one can speculate on models in which 
significant inhomogeneities in r], (Ary/ry^l), may occur. Such inhomogeneities, for example, 
might arise in nonadiabatic (thin wall) models [9] whenever the velocity of an expanding 
broken phase bubble varies during the transition. In nonadiabatic scenarios the rate of 
baryogenesis depends strongly upon the velocity of the wall. We note that this velocity 
may change at only about the 10% level during the course of the transition. However, 
there remains considerable uncertainty in the determination of the wall velocity in these 
models. As recently pointed out, this effect may occur in adiabatic models as well [10]. 
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Another possibility might be the formation [11] of distinct domains of baryon- number 
and anti-baryon number. After annihilating they could leave behind a small number of 
baryon bubbles containing all of the net baryon number. Finally, strong spatial inhomo- 
geneities may result from any scenario in which most of the generated net baryon number 
is associated with the collisions of bubble walls at the end of the transition. Although such 
models are speculative it is nevertheless interesting to investigate the constraints which 
might be placed on such scenarios from primordial nucleosynthesis. 

If fluctuations occur, the bubble size at coalescence will probably provide a typical 
length scale of fluctuations in rj, but fluctuations can occur on larger scales than that. 
However, it is difficult to quantify that length scale. Thermal and/or quantum nucleation 
is especially difficult to follow at the electroweak epoch because the nucleating action may 
be dynamically renormalized by the presence of bubbles of broken phase [12]. Another 
complication may be hydrodynamic instability of phase boundaries [13]. 

Despite these caveats it is nevertheless instructive to consider simple models of ho- 
mogeneous nucleation of phase in the small supercooling limit [14] . In these models the 
nucleation rate per unit volume is assumed to be, 

p(T) ^ CT^e-^^^) , (1) 

where S{T) = a{T)(Tc/{Tc — T)) is the nucleating action, with a{T) a monotonically 
increasing function of temperature, and where C is a scale factor of order unity. Integrating 
the nucleation rate through the end of the phase transition, and assuming that bubble walls 
move at the speed of light, yields an estimate for the time required for bubbles to coalesce. 



We can express this coalescence time (or bubble size at coalescence) as a fraction 5 of the 
Hubble time (or horizon scale) H~^ [13], 

5^(4Sln(^)) \ (2) 

where B is the logarithmic derivative of the nucleating action S, in units oi H~^ at the 
epoch of the phase transition. The value of B depends on calculable details of models 
for the electroweak transition and is within one or two orders of magnitude of unity. The 
horizon size is, 

^ UTrsJ ^ T2 ^^^-^^'^'^HlOO; llOOGeVy' ' ^^^ 

where g = X^bfi'b + '^/^'l2f9f is the total statistical weight in relativistic bosons (^r^) 
and fermions (gf) at temperature T, and nipi is the Planck mass. In the standard model, 
g ^ 100 for an electroweak transition at T = 100 GeV. The total statistical weight is slightly 
uncertain due to the unknown top quark mass and extra degrees of freedom associated with 
extensions of the standard model. 

The average bubble size at coalescence 6H~^ is a result of competition between the nu- 
cleation rate and the very slow expansion of the universe. Most bubbles will have size 5H~^ 
at coalescence. This follows on noting that larger bubbles would have to be nucleated early, 
near Tc, where the nucleation rate is exponentially suppressed. Smaller bubbles would have 
to be nucleated near the end of the phase transition where the effective nucleation rate 
is again small, since very little unbroken phase would remain. A nucleation/coalescence 
epoch which approximates the homogeneous nucleation scenario will leave a nearly regular 
lattice of bubbles at coalescence. 



Even though several electroweak scenarios have been proposed [8,9,15] not much is 
known about the actual nucleation scale SH~^ and the expected coexistence temperature 
Tc in these models. It has been argued that the minimal standard model gives 6 ~ 
10~^ [16]. These models have not been investigated in sufficient detail to ascertain the 
relationship between the nucleation scale 5H~^ and the scale of separation between centers 
of fluctuations in r], which we shall denote 5t^H~^. However, it is possible that 5^ ^ 5, 
corresponding to less than, or equal to, one fluctuation produced per nucleated bubble. 

In Ref.[17], two of us (hereafter referred as JF) have studied in detail the evolution 
of fluctuations from T ^ 100 MeV to T ^ 1 keV taking into account neutrino, photon, 
and baryon dissipation processes. JF's results indicate that fluctuations generated at 
the electroweak epoch may survive through the epoch of primordial nucleosynthesis. If 
fluctuations with particular characteristics produced at an early epoch did survive, their 
presence could alter the nuclear abundance yields emerging from primordial nucleosynthesis 
[18,19]. If these abundance yields do not agree with observationally inferred primordial 
abundances, then we can conclude that these fluctuations could not have existed. This, in 
turn, would allow us to constrain the fluctuation generation mechanism. 

We follow JF and deflne the amplitude of fluctuations A(x), in terms of the spatial 
distribution of baryon-to- photon number, r]{x), and its horizon average, fy, by r]{x) = f]{l + 
A(a;)). The corresponding distribution in entropy-per-baryon is then s(a;) = s(l+A(x)) , 
where the average conserved entropy-per-baryon in units of Boltzmann's constant is s ~ 
2.63 X lO^O^ h~'^ . In this expression O^ is the fraction of the closure density contributed 



by baryons and h is the present Hubble parameter in units of 100 km s~^ Mpc~^. In this 
paper A{x) always refers to the initial amplitude of the fluctuations. 

The primary criterion for fluctuation survival is that the scale associated with the 
separations of the centers of fluctuations (5^if~^) be comparable to, or exceed, the comov- 
ing proton diffusion length (dioo) at the beginning of the nucleosynthesis epoch [17,18,19]. 
Here dioo is the comoving proton diffusion length referenced to the epoch of T = 100 GeV 
(see for example eq.4). Were this condition not satisfled, baryon diffusion would erase 
fluctuations in r] prior to nucleosynthesis. The proton diffusion length is actually a fairly 
sensitive function of amplitude (1 + A). In Figure 1 we give the comoving proton diffusion 
length dioo at the epoch T = 500 keV as a function of (1 + A). This temperature very 
roughly corresponds to the epoch of weak freeze-out, where the neutron-to-proton intercon- 
version rate from lepton capture falls below the free neutron decay rate. Note that higher 
baryon density implies a smaller diffusion length for baryons. Whenever (1 + A)^10^ the 
baryon diffusion length corresponds to dioo ~ 0.1 cm. 

We also can describe fluctuations by their separation length scale, /loo, where we 
express a length scale co-moving with the Hubble expansion in terms of its proper length 
at an epoch where T = 100 GeV. The corresponding proper length at any epoch where 
the temperature is T is then 

where R and -Rioo stjce the scale factors at an epoch of temperature T and 100 GeV, 
respectively, Tioo = 100 GeV, and where g and gioo are the statistical weights in relativistic 



particles at an epoch of temperature T and 100 GeV, respectively. In this expression we 
have assumed that the co- moving entropy density is conserved. 

In order for a fluctuation to affect the outcome of nucleosynthesis /loo ^ ^locT ~ '^loo- 
This scale is found from detailed nucleosynthesis calculations to be roughly the scale of the 
proton diffusion length at the nucleosynthesis epoch. Physically, the origin of this limiting 
length is that any fluctuation scale smaller than the proton diffusion length will be damped 
out by baryon diffusion prior to nucleosynthesis. Therefore, the minimum fluctuation scale 
for inhomogeneous nucleosynthesis effects can be expressed in terms of a fraction of the 
horizon scale H~^ at any epoch as 

r^^ f8Tr^\ 1/3 1/6 ^^100 ., X 

drmn - ^_, ~ ^100 ^ qq J ^100^ ^^^ , (5«J 

5_ . (6.9 X 10-) r^) '^' (^) ''' (^^) (^^) . (56) 

^ ^VlOO/ \gim^ V 0.1cm/ VlOOGeV/ ^ ' 

Note that 5min ^ 1 for T ^ 1.45 TeV. We conclude that micro-physical, subhorizon-scale 
fluctuation-generating processes operating at epochs for which T ^ 1.45 TeV conceivably 
could have constrainable nucleosynthesis signatures. Fluctuations in tj on initially super- 
horizon scales / which satisfy / ^ Imin are similarly at risk of running afoul of primordial 
abundance constraints. We note, however, that even fluctuations with I "zlmin may yet 
survive to affect nucleosynthesis if they have amplitudes large enough that the length scales 
of their high density regions, /f^g? exceed (iioo- In this case baryons would be unable to 
diffuse out of the high density cores of fluctuations prior to nucleosynthesis. 

In any scheme for baryogenesis associated with an electroweak symmetry breaking 
epoch at temperature T we must produce the average proper baryon number density 



within the horizon 



s 
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where the entropy per unit proper volume is 5" ~ {2n'^/45)gT^. Homogeneous and in- 
homogeneous standard big bang nucleosynthesis calculations together with observational 
abundance constraints imply that O5 ^ 0.01h~^ [1,18]. 

Assume that baryons are distributed in high density regions with baryon number 
density n^, which in total occupy a fraction fv of the horizon volume, and in low density 
regions with baryon number density n^. In this case, we can write 

Ub = fv< + (1 - fv)nt (7) 

We define A^ = n^ /fib and A^ = n^ /fib, so that the density contrast between high 
and low density regions is A = Ah/A^. If the horizon is filled with a regular lattice of 
fiuctuation cells whose centers are separated by /fgo, then the length scale of high density 
regions is /foo = fy ^foo- 

In Ref.[18], three of us have calculated in detail the outcome of primordial nucleosyn- 
thesis with inhomogeneous initial conditions. In these calculations the nuclear reaction 
rates were coupled to all significant fiuctuation dissipation processes: neutrino heat trans- 
port, baryon diffusion, photon diffusive heat transport, and hydrodynamic expansion with 
photon-electron Thomson drag. The light element abundance yields are found to be in- 
consistent with observations for all but a very narrow range of fiuctuation characteristics. 
This is why nucleosynthesis is so powerful in constraining primordial inhomogeneities. 
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A representative case of these calculations is displayed in Figure 2. In this figure we 
show the ^He mass fraction and number fraction of ^Li emerging from an inhomogeneous 
big bang with fluctuation separations /foo- Dotted lines indicate lYoo' (BDL) and the elec- 
troweak horizon scale (EWH). Though we show results for particular values of fluctuation 
amplitude and gaussian width aioo (roughly, fy ~ ('^loo/^foo)"^)' ^^^ flgure illustrates some 
general trends for abundance yields as a function of length scale. In particular, we note 
that when /^q* ;$ IIqq ^ EWH the abundances of ^He and/or ^Li always exceed observa- 
tional limits [20]. This is a general feature of inhomogeneous primordial nucleosynthesis 
whenever I'Iqq is below the electroweak horizon scale. It is, however, intriguing that the 
electroweak horizon is close to the minimum in ^He and ^Li (the "helium dip"). 

We have also explored inhomogeneous nucleosynthesis yields for the light elements as 
functions of initial density contrast and volume flUing fraction, A and fy, respectively. 
Figure 3 shows the results of numerous numerical calculations for I'Iqq = 0.5 cm. In this 
flgure the parameter space of A and fv laying to the right of the shaded line gives ^He 
overproduction (^He mass fraction > 24%). In general, we flnd that ^He and/or ^Li are 
overproduced whenever 10% or more of the baryons reside in high density regions and 
either /foo > lYoo' or I^^qq > <iioo- 

This constraint can be put in the context of electroweak baryogenesis with a simple 
model. Assume that an electroweak phase transition has produced a regular lattice of 
bubbles at coalescence, all of equal size. In fact, we expect a distribution of bubble sizes, 
but for now assume all bubbles have size equal to the nucleation scale 5H~^. Since we 
expect 5 ~ 10""^ there will be roughly 5~^ ~ 10^ bubbles within the horizon. As they 
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expand toward coalescence, these 10^ bubbles must produce the horizon-averaged baryon 
number density fib- It is likely, however, as noted above that not all bubbles will contribute 
equally to this average. Assume that a fraction / of the bubbles contributes a substantial 
fraction b of the total baryon number in the horizon. As far as the net baryon distribution 
is concerned, this scenario approximates a regular lattice of fluctuations with effective 
horizon-fraction separations 5,, ~ f~^'^5. 

Equation (7) shows that the baryon distribution is characterized by two independent 
quantities, which we take to be fv and A. In the above hypothetical scenario it is clear 
that fv < /, with equality obtaining when the baryon number produced in a bubble is 
uniformly distributed over the volume swept out by the bubble wall. The density contrast 
will be 

assuming that each of the fd~^ bubbles contributes equally to bfib- Equality in equation 
(8) obtains in the limit of uniform baryon distribution across each of the "significant" /<^~^ 
bubbles. Note that we also assume that significant bubbles are uniformly distributed in 
space. In this model, b = (/yn^)/nb = fv^H- Since Ah = 1 -|- A, we may conclude that 
(1 + A) ^ b/ fv- The total fraction of baryons in high density regions is /vA, which is just 
b in the limit where A/y << 1. Our nucleosynthesis constraint will apply whenever more 
than 10% of the baryons are in high density regions, b ^ 0.1, and either the fiuctuation 
separation exceeds the minimum value required for nucleosynthesis effects, f~^^^6 > Smim 
or the high density region length scale exceeds the proton diffusion length, I^qq ^ c^ioo- 
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For example, in the context of this model assume that only 1000 bubbles out of the 
total of 10^ bubbles in the horizon produce 90% of the baryon number. This implies 

-1 /o 

that /y Ri / fti 10~^ and dr, ^ fy S ^ 10^5 ^ 0.1 > dmin- This last comparison 
follows on assuming that S ~ 10~^. Equation (8) implies that A si 9 x 10^. Since 
b = 0.9 exceeds 0.1 and 5^ > dmin we would conclude that this scenario is incompatible 
with nucleosynthesis constraints and, therefore, ruled out. In general, whenever / < 
(S/dmin)^ ~ 3 X 10~^(100GeV/T)^ and b > 0.1 the nucleosynthesis constraint wiU be 
violated. In other words, any scenarios where fewer than about 3000 bubbles contribute 
more than 10% of the baryon number can be ruled out. 

In some models of electroweak baryogenesis the significant bubbles for baryon pro- 
duction might be those which are nucleated earliest. We would then expect these bubbles 
to be larger at coalescence than the average bubble {6H~^). If the ratio of significant 
bubble size to the average bubble size is r, and there are fS~^ significant bubbles, then the 
effective volume filling fraction for the higher density regions of the baryon distribution is 
roughly fv ^ r^f- Note that in this case, however, the effective horizon-fraction separation 
is still 5^ ^ f-^/^S. 

Though the detailed relationship between the fluctuation scale and amplitude and 
such model parameters as the Higgs particle mass, the top quark mass, and the critical 
temperature T^ are as yet poorly understood, they are in principle calculable. The param- 
eters of electroweak baryogenesis models must not lead to a violation of the constraints on 
fluctuation characteristics as describes above. This is probably readily attainable for some 
models, but may represent a restriction for others. 

12 



Acknowledgements 

We wish to thank N. J. Snyderman and H. Kurki-Suonio and M.S. Turner for many 
useful conversations. We also wish to thank an anonymous referee for many useful com- 
ments regarding a preliminary version of this manuscript. This work was supported in part 
by NSF Grant PHY91-21623, IGPP LLNL Grant 93-22, and a California Space Institute 
Grant. It was also performed in part under the auspices of the U.S. Department of Energy 
by the Lawrence Livermore National Laboratory under contract number W-7405-ENG-48 
and DoE Nuclear Theory Grant SF-ENG-48 and by the University of Chicago under DOE 
grant FG02-91ER40606. 



13 



References 

[1] E. W. Kolb and M. S. Turner, " The Early Universe " , Addison- Wesley (1990). 

[2] A. D. Sakharov, JETP Lett. 5 (1967) 24. 

[3] V. A. Kuzmin, V. A. Rubakov, and M. E. Shaposhnikov, Phys. Lett. 155B (1985) 36; 

E. W. Kolb and M. S. Turner Mod. Phys. Lett. A2 (1987) 285. 

[4] M. E. Shaposhnikov, Phys. Lett. B277 (1992) 324. 

[5] M. Dine, P. Huet, R. Singleton, and L. Susskind, Phys. Lett. 257B (1991) 351; M. 

Dine, O. Lechtenfeld, B. Sakita, W. Fischler, and J. Polchinski, Nucl. Phys. B342 (1990) 

381. 

[6] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, Phys. Lett. 245B (1990) 561; A. G. 

Cohen, D. B. Kaplan, and A. E. Nelson, Phys. Lett. B263 (1991) 86; A. G. Cohen, D. B. 

Kaplan, and A. E. Nelson, Nucl. Phys. B349 (1991) 727. 

[7] L. McLerran, M. Shaposhnikov, N. Turok, and M. Voloshin, Phys. Lett. B256 (1991) 

451; N. Turok and J. Zadrozny, Nucl. Phys. B369 (1992) 729. 

[8] A. G. Cohen, D. B. Kaplan, and A. E. Nelson, Ann. Rev. Nucl. and Part. Sci. in 

press. 

[9] A. E. Nelson, D. B. Kaplan, and A. G. Cohen, Nucl. Phys. B373 (1992) 453. 

[10] M. Joyce, T. Prokopec, and N. Turok, preprint |hep-ph/9401352 . 



[11] D. Comelh and M. Pietroni, Phys. Lett. B 306 (1993) 67; J.R. Espinosa and M. 
Quiros, Phys. Lett. B 307 (1993) 106; D. Comelh, M. Pietroni, and A. Riotto, Nucl. 
Phys. B 412 (1994) 441. 

[12] M. Gleiser and E. W. Kolb, Phys. Rev. Lett. 69 (1992) 1304. 

14 



[13] K. Freese and M. Kamionkowski, Phys. Rev. Lett. 69 (1992) 2743; P. Huet, K. 

Kajantie, R. Leigh, B. Liu, and L. McLerran, preprint SLAC-PUB-5943. 

[14] C. Hogan, Phys. Lett. 133B (1983) 172; S. Coleman, Phys. Rev. D15 (1977) 2929. 

[15] M. Dine, R. G. Leigh, P. Huet, A. Linde, and D. Linde, Phys. Rev. D46 (1992) 550; 

K. Enqvist, J. Ignatius, K. Kajantie, and K. Rummukainen, Phys. Rev. D45 (1992) 3415; 

B. H. Liu, L. McLerran, and N. Turok, Phys. Rev. D46 (1992) 2668. 

[16] G. W. Anderson and L. J. Hall, Phys. Rev. D45 (1992) 2685. 

[17] K. Jedamzik and G. M. Fuller, Astrophys. J. 423 (1994) 33. 

[18] K. Jedamzik, G. M. Fuller, and G. J. Mathews, Astrophys. J. 423 (1994) 50. 

[19] J. H. Applegate, C. J. Hogan, and R. J. Scherrer, Phys. Rev. D35 (1987) 1151; C. 

R. Alcock, G. M. Fuller, and G. J. Mathews, Astrophys. J. 320 (1987) 439; See the review 

by R. A. Malaney and G. J. Mathews, Physics Reports, 229 (1993) 145. 

[20] The mass fraction of ^He, Yp, must satisfy Yp ^0.24. It is believed that the primordial 

^Li abundance is the Population II result ^Li/Hsi lO"^'^, but it is certainly smaller than 

the Population I value of ~ 10~^. See discussions in A. M. Boesgaard and G. Steigman, 

Ann. Rev. Astron. Astrophys. 23 (1985) 319; Walker et al., Astrophys. J., 376 (1991) 51; 

and ref. [18]. 



15 



Figure Captions 

Figure 1: The co-moving proton diffusion lengtli, dioo, at temperature T — 500 IceV 
as a function of fluctuation amplitude (1 + A). Ttie calculation assumes O^/i^ = 0.0125. 

Figure 2: The ^He mass fraction (upper panel) and ^Li number fraction relative 
to hydrogen from inhomogeneous nucleosynthesis calculations are shown as functions of 
the CO- moving fluctuation separation /fgo- These calculations employ an initial density 
contrast A = 1.25 x 10^, Gaussian width aioo = ^foo/^O, and O^/i^ = 0.0125. The dotted 
lines indicate the minimum scale above which fluctuations survive baryon diffusion, BDL, 
and the electroweak horizon scale, EWH. 

Figure 3: Constraints from inhomogeneous big bang nucleosynthesis on models with 
a high-density region volume fllling fraction fy and density contrast A = n^ /n^. Any 
models with parameters falling to the right of the shaded line will overproduce ^He relative 
to observational limits. For this calculation we assumed spherical square- wave fluctuations 
with fluctuation mean separation /Iqq = 0.5 cm and O^/i^ = 0.0125. 
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